arXiv:1506.06292v2 [math-ph] 25 Aug 2015 


Decoherence of a qubit as a diffusion 
on the Bloch sphere 

Katarzyna Siudzinska* and Dariusz Chruscinski 
Institute of Physics, Nicolaus Copernicus University, 
ul. Grudzi§dzka 5, 87-100 Toruh, Poland 

August 26, 2015 


Abstract 

We analyze qubit decoherence in the framework of geometric quantum mechanics. In this 
framework the qubit density operators are represented by probability distributions which 
are also the Kahler functions on the Bloch sphere. Interestingly, the complete positivity of 
the quantum evolution is recovered as ellipticity of the second order differential operator 
(deformed Laplacian) which governs the evolution of the probability distribution. 


1 Introduction 


In the geometric approach to quantum mechanics quantum states are represented not 

by vectors in the Hilbert space 77, but as points in the projective Hilbert space P77. A point in 
P77 corresponds to a ray in 77 passing through U or, equivalently, to the rank-1 projector |U)(UI 
(see also mum and [8] for recent reviews). If 77 = then 

P77 = CP" = U{n + 1)/U(n) . (1) 

It turns out that CP" defines an n-dimensional complex space equipped with rich geometrical 
structures: the Fubini-Study metric g and the symplectic form w, such that 


}C = g + iuj. 


( 2 ) 


The triple (CP",g',a;) defines a Kahler space (see the recent paper [9]). Now, for any operator 
A G H(77) one defines a function Ja '■ CP" —)• C via 
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fAim 


(U|A|U) 

(UIU) 


(3) 


1 




These functions are called the Kahler functions [lOj and form a linear subspace is the space of all 
functions T{CP^) := {/ : —> C}. An equivalent definition of the Kahler function uses the 

following property: for any function / : CP” —>■ C one dehnes the corresponding Hamiltonian 
vector field Xf hy d f = u;(Aj, •). Now, / is Kahlerian iS Xf is a Killing vector field, that is, 
^Xfd = 0. Note that if A,P G B{T-i) and C = —2i[A,P], then for the corresponding Kahler 
functions fA, fB, fc one has 

fc = {fA,fB} = -XfMB) , (4) 

where {/a,/b} := denotes the Poisson bracket on CP” given by the symplectic 

form uj. The corresponding Hamiltonian vector fields satisfy Xf^ = [Xf^,Xfg]. Finally, the 
adjoint map adn ■ B{'H) -A B{'H), 


aduiA) =-\[H,A] , ( 5 ) 

gives rise to the map P(CP”) i-)- P(CP”) defined by a ^ 4{/i, a} = —^Xh{a). 

For the qubit case the corresponding complex projective space CP^ is nothing but the Bloch 
sphere. Introducing € C^ by 

Wo = cos - , = sm-e 

with 9 G [0, tt) and 4> G [0, 27r) one finds the Fubini-Study metric 

ds^ = ^(d0^ + sin^ 0d(^^) , (6) 

and the symplectic 2-form 

w = ^ sin0d0 A d(^ . (7) 

It defines a volume element on CP^ (c.f. i) 

Vol(Cpi) = [ u} = 7r . (8) 

Jcp^ 

Taking the basis in M 2 (C) of {(To = I,ui,(T 2 ,us}, one finds for fk{['<P]) = {'<PWk\'4’)' 

/o( 6 »,</)) = !, /i( 6 »,^) = sin0cos(?!) , f2{9,(j)) = sm9sm(j), / 3 ( 6 », ^) = cos6» . (9) 

These functions provide a basis in the space of Kahler functions on CP^. The corresponding 
vector fields are Xq = 0 and 


Xi = —4 sin (pdg — 4 cot 9 cos 

A 2 = 4cos (/)5e — 4cot 0 sin(^(90, (10) 

^3 = 49^ , 
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Note that reproduce three components of angular momentum, and hence 

+ Lg = , (11) 

where A is the Laplacian on CP^. The Kahler functions {fi, f 2 i fa} a-re dipole spherical har¬ 
monics, therefore Afk = —8/^0 

In this paper we consider the above geometric structures to provide a description of deco¬ 
herence processes on CP^. Section [2] provides a simple model of qubit decoherence, which is 
then analyzed within the geometric approach in Section [3l Section [3] discusses the conditions 
for complete positivity, and the final conclusions are collected in the last section. 


2 Qubit decoherence 


The Markovian master equation has the following form mm, 


Pt = - i[H, Pt] + - ^ Ik 


VkPtVl--[VlVk,pt]+ 


( 12 ) 


where H denotes an effective Hamiltonian of the system, 14 are the noise operators, and the 
coefficients 7 fc > 0 {[A,B]^ = AB + BA). Now, consider the following master equation for the 
evolution of a qubit. 



3 

^ ^ ^kipkPt^k Pt) ) 
k=l 


(13) 


where 7 fc > 0 denote the decoherence rates, and ak are the Pauli matrices. The most convenient 
way to represent a solution of (1131) is to use the Bloch representation. 


Pt = ]^\^ + '^Xk{t)akj ■ (14) 

One easily finds that the Bloch vector x(t) = (xi(t),X 2 (t),X 3 (t)) satisfies 

xi{t) = -[72 + lalxiit) , X 2 {t) = -[71 -h 'y3]x2{t) , X3{t) = -[71 -h 'y2]x3{t) , (15) 

and hence 

xi{t) = , X 2 {t) = , X 3 {t) = e“['’'^’'''^^]*X 3 ( 0 ) . (16) 


The above solution represents an anisotropic decoherence of a qubit. Clearly, in the isotropic 
case, that is, if 71 = 72 = 73 = 7, one has x(t) = e“^'^*x(0), which describes a depolarizing chan¬ 
nel. Asymptotically Xk{t) —)• 0, which means that pt approaches a completely decohered 

maximally mixed state. 

^Dipole functions Yi^ {m = —1,0,1) for the Laplacian on the unit sphere satisfy AYim = —2Yim. Here, due 
to the factor in da one has ‘8’ instead of ‘2’. 
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The evolution p ^ pt may be described by the following completely positive trace-preserving 
map, 


= '^Pa{t)(7apaa 


(17) 


Q :=0 


with the probability distributions Pa{t) defined by |14j : 


Plit) = 

b 

_ g-[7l+72]t _|_ g-[72-l-73]t _ g-[73-l-7l]i 

P2{t) = 

b 

_ g-[7l+72]i _ g-[72-l-73]t _j_ g-[73-l-7l]i 

Psit) = 

b 

_|_ g-[7l+72]t _ g-[72-l-73]t _ g-[73-l-7l]i 


and po{t) = 1 — Pi(t) — P 2 {t) — Psit)- Note that the Bloch vector x(t) stays within the Bloch 
ball whenever 


7i + 72 > 0 , 72 + 73 > 0 , 73 + 7i > 0 , (18) 

which is a mnch weaker requirement than 7i, 72 , 73 > 0 responsible for complete positivity of 
quantnm evolution. Positive decoherence rates correspond to the Markovian evolntion, whereas 
the decoherence rates satisfying (|18p describe the non-Markovian evolntion M- We shall discnss 
it in more details in Section 01 

Observe that (jl.lh may be equivalently rewritten as follows, 


Pt = - 


1 

4 


3 

'^lk[o-k,Wk,Pt]] ■ 
k=l 


(19) 


Dne to the “double commutator” structure this equation is perfectly suited for the geometric 
reformulation and it provides the starting point of the geometric analysis. 


3 Geometric description 

For a given density operator p let us define 


P([V’]) := 


2 (V>|p|-0) 

TT (V'1'0) 


It is clear that p(['i/’]) is a probability distribntion on CP^, that is, p(['0]) P 0 and 


( 20 ) 



( 21 ) 


where uj is defined in d?]). It corresponds to a legitimate density operators if and only if p is a 
Kahler function. Note that p describes a pure state p = |V’)(V’I iff p{W\) — f • 
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For any unitary operator U : T-L ^ T-L one defines a mapping U : J^(CP"') as 

follows, 

= /([^'V]) • (22) 

If pt satisfies the von Neumann equation 

Pt = - i[H, Pt] , (23) 

then the corresponding family of probability distributions Pt([V’]) satisfies 

Pt = ^{h,pt} =-Xh Pt , (24) 

and the solution is given by 

Pt([V’])=Po([e*^V]) • (25) 

This solution corresponds to the rigid rotation on CP^. For H = uqI + the initial 

probability distribution po is rotated along n by an angle |n|t. Now, consider the dissipative 
evolution of pt given by ([T9]). It induces the following dynamical equation for pt, 

1 ^ 

Pt = — X] 5'fc^fcPt • (26) 

° k=l 

Interestingly, in the isotropic case 7i = 72 = 73 = 7 one finds 

3 

Y,lkXl = A^X, (27) 

k=l 

where A denotes the Laplacian on CP^, and hence (I26p reduces to the diffusion equation on the 
Bloch sphere, 

Pt = '^Xpt . (28) 

The diffusion equation on a sphere (I28p may be easily solved by expansion of the probability 
distributions into harmonic functions Yim- Assuming that 


00 I 


PO — ^ ^ ^ ^ tllmYlm 

(29) 

1=0 m=—l 


and taking into account 


XYim = —41(1 + l)l)m ) 

(30) 

one hnds the solution in the following form, 


00 1 

Pt = ^ ^ . 

1=0 rn=—l 

(31) 


5 


Asymptotically pt —)■ aoo = Now, if pt represents a legitimate quantum state, then the 
monopole and dipole functions are the only functions to enter the expansion (I29p . and therefore 

Pt = - (l + + X 2/2 + X3/3]) , (32) 

TT 

where x = {xi^X 2 ,x^) denotes the Bloch vector with the initial state pQ. Note that the Bloch 
vector X corresponding to the pure state po = IV’o)(V'ol satisfies |x| = 1, and hence may be 
represented by x = (sin0ocos0o,sin0osin(?!)o,cos0o) for some point {9q,4>q). Moreover, 
implies 

Pt(0o>o) =-(1 + e-"^*) , (33) 

TT 

which shows that pt{do, </>o) evolves from ^ at t = 0 to ^ as t 00 . Clearly, pt{do, cpo) > pt{9, </>) 
for any {9,(j)) / (^ct/'o)- 

In the anisotropic case finding the solution of (1261) is much more difficult since are no 
longer the eigenfunctions of Ylk=i ^ > 1- However, for the Kahler functions pt one finds 

Pt = \{^ + , (34) 

which, of course, reproduces (fTHl) . In the axial case, i.e. 71 = 72 = 0 and 73 = 7 , one has 

Pti9,(j)) = jdlpt{9,(j)) , (35) 

which is the diffusion equation on the 1 -dimensional torus parameterized by (p and corresponds 
to the phase damping channel. 

4 Positivity vs. complete positivity 

Let us recall that the positivity of jk provides a necessary and sufficient condition for complete 
positivity of the solution p ^ pt = Af[p] of the original master equation (fT^ . However, this 
condition is somehow missing when one considers the evolution of pt on the Bloch sphere. Note 
that in order to have a legitimate probability distribution provided by the Kahler function pt 
one needs to impose weaker conditions defined in (jl 8 l) . Actually, these conditions guarantee 
that p 1 -^ Pt = At [p] is a positive map for all t > 0 [H] . 

It should be stressed that equation ([M]) . when supplemented by (fTHl) . guarantees that pt is a 
probability distribution on CP^ if po is a Kahler function. However, if pQ is not Kahlerian (which 
means that it contains the multipole expansions of harmonics higher than dipole functions), then 
Pt needs not to be a probability distribution for all t > 0. 

Let us introduce A..y := 4 7 fcAi| (one may call it an anisotropically deformed Laplacian). 

Theorem 1 The operator is elliptic iff 'jk ^0 for k = 1,2, 3. 
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Proof: The anisotropically deformed Laplacian can be rewritten as 


=4 ^ 7 i sin^ (f) + "i 2 cos^ (j^d^ + A ^71 cot^ 9 cos^ </> + 72 cot^ 6 sin^ cj) + 73 ^ d\ 
+ 4 cot 9 sin(2(/)) ^71 — 72 ^ ded^f) + 4 cot 9 ^71 cos^ 0 + 72 sin^ do 

[^72 - 7i)'9<i.- 


sin(2(;i)) 


sin 


3 + cos 


Now, is elliptic on CP^ iff the following diffusion matrix 


with the coefficients 


A{9A) 


aii { 9 , 4 >) 012(0,0) 

012(0,0) 022(0, 0 )J ’ 


( 36 ) 


(37) 


aii(0, 0) = 4^71 sin^ 0 + 72 cos^ 0^ , 

0 : 22 ( 0 , 0) = 4^71 cot^ 0 cos^ 0 + 72 cot^ 0 sin^ 0 + 73 ^ , (38) 

ai 2 ( 0 , 0 ) = 2 cot 0 sin( 20 )(^ 7 i - 72 ^ , 

is positive definite for each 0, 0. Positivity of A{9, 0) is equivalent to the conditions 

oii(0,0)>O, 022 ( 0 , 0 ) >0, 011 ( 0 , 0 ) 022 ( 0 , 0 ) > 012 ( 0 , 0 ) . (39) 

Now, oii(0,0) > 0 iff 71 > 0 and 72 > 0, and the condition oii(0, 0 ) 022 ( 0 , 0) > 0 ^ 2 ( 0 , 0 ) is 
equivalent to 


7172 cot^ 0 + 7173 sin^ 0 + 7273 cos^ 0 > 0 , (40) 

which finally implies 73 > 0. □ 

Corollary 1 Equation i26\) provides a legitimate description of the pt corresponding to a com¬ 
pletely positive evolution of pt if and only if Ay is an elliptic operator which means that 
defines the Fokker-Planck diffusion equation on CP^. 

5 Conclusions 

We provide the description of qubit decoherence within geometric approach to quantum me¬ 
chanics. It is shown that Markovian master equation for density operator is replaced by the 
diffusion-like equation for the probability distribution on the Bloch sphere. It turns out that 
complete positivity of the evolution is equivalent to the requirement that the evolution of the 
probability distribution is governed by the legitimate Fokker-Planck equation m- In the forth¬ 
coming paper we generalize the presented results to the decoherence processes on CP”. 
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The above analysis may be easily generalized to the time-dependent case, i.e. when 7fc(t) do 
depend on time. In this case, the diffusion matrix A[9, (/>; t) defined in (1371) is time-dependent 
(via the time-dependence of 7 fc(t)). One may characterize the Markovianity of the evolution by 
analyzing the properties of A(6,(p;t). The non-Markovian quantum evolution attracts recently 
considerable attention (see [MHZ] for the recent reviews). It turns out that memory effects 
might play an important role in modern quantum technologies and quantum information. Now, 
the evolution is CP-divisible |18] if and only if the matrix A{6, (/>; t) is positive-definite for all 
t > 0 and it is P-divisible if and only if A{6, (p; t) is positive-definite on the subspace of Kahler 
functions. Hence, in the geometric approach Markovianity of the evolution is controlled by the 
time-dependent diffusion matrix defining the Fokker-Planck equation. 
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